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We investigate the transverse dynamics in Au+Au collisions at ^snn = 200 GeV by emphasis 
upon the interplay between soft and hard components through pt dependences of particle spectra, 
ratios of yields, suppression factors, and elliptic flow for identified hadrons. From hydrodynamics 
combined with traversing minijets which go through jet quenching in the hot medium, we calculate 
interactions of hard jets with the soft hydrodynamic components. It is shown by the explicit 
dynamical calculations that the hydrodynamic radial flow and the jet quenching of hard jets are 
the keys to understand the differences among the hadron spectra for pions, kaons, and protons. 
This leads to the natural interpretation for Np/N^ ~ 1, RaaS^ 1 for protons, and > 1)% recently 
observed in the intermediate transverse momentum region at Relativistic Heavy Ion Collider (RHIC). 



PACS numbers: 24.85.+p,25.75.-q, 24.10.Nz 



I. INTRODUCTION 

A vast body of data has already been collected and 
analyzed during past few years at Relativistic Heavy Ion 
Collider (RHIC) Q toward a complete understanding of 
the dense QCD matter which is created in high energy 
heavy-ion collisions. 

At collider experiments, it is well known that high 
perturbative QCD (pQCD) processes become so large as 
to observe jet spectra. One of the most important new 
physics revealed in heavy ion collisions at RHIC ener- 
gies is to study propagation of (mini-)jets in dense QCD 
matter. Jet quenching has been proposed as a possi- 
ble signal of deconfined nuclear matter, the quark gluon 
plasma (QGP) (for a recent review, see Ref. 0). Over 
the past years, a lot of work has been devoted to study 
the propagation of jets through QCD matter 0, 0,0,01 . 

Recent data at RHIC indicate that both the neutral 
pion [U, and the charged hadron 0, 0, spectra 
at high pt in central Au+Au collisions are suppressed rel- 
ative to the scaled pp or large centrality spectra by the 
number of binary collisions. However, protons do not 
seem to be quenched in the moderate pt range fl3j . Fur- 
thermore, the proton yield exceeds the pion yield around 
p T 2-3 GeV/c which is not seen in elementary hadronic 
collisions [T^]. The STAR Collaboration also shows that 
A/K° ~ 1 at a transverse momentum of 2-3 GeV/c [l4j . 
pQCD calculations are successful in describing hadron 
spectra in Au+Au collisions as well as pp collisions by 
taking account of nuclear effects such as Cronin effect, 
nuclear shadowing effect, and energy loss of jets [l5j . 
However, large uncertainly of the proton fragmentation 
function makes the understanding of the baryon produc- 
tion mechanism unclear [T^ | even in pp collisions. On the 
other hand, several models have been proposed by con- 
sidering interplay between non-perturbative soft physics 
and pQCD hard physics: baryon junction [TtI Ifs|. par- 



ton coalescence 0, |2(], El I22II23I , medium modification 
of the string fragmentation |24j . and a parametrization 
with hydrodynamic component combined with the non- 
thermal components in order to explain the anoma- 
lous baryon productions and/or large elliptic flow discov- 
ered at RHIC. 

It is said that hydrodynamics j2|| [27|, l28l |29| works 
very well for explanation of elliptic flow data at RHIC 
energies, in the low p^ region, in small centrality events, 
and at midrapidity, including the mass dependence of 
hadrons (for recent reviews, see Ref. This suggests 

that hydrodynamics could be reliable for the description 
of the time evolution of soft sector of matter produced 
in high energy heavy ion collisions at RHIC. Certainly, 
it is more desirable to describe the time evolution of the 
whole stage in high energy heavy ion collisions by simu- 
lating collisions of initial nuclear wave functions. Instead, 
they simply assume that the system created in heavy ion 
collisions reaches local thermal equilibrium state at some 
time. 

Due to the above two reasons, a model which treats a 
soft sector by hydrodynamics and a hard sector based on 
a pQCD parton model is turned out to be useful in or- 
der to understand experimental data at RHIC from low 
to high px- Indeed, first attempts based on this concept 
has been done by pQCD calculations which include hy- 
drodynamic features 0, ■ Motivated by these 
works, we have recently developed a two component dy- 
namical model (hydro+jet model) l34l with a fully three 
dimensional hydrodynamic model [28| for the soft sector 
and pQCD jets for the hard sector which are computed 
via the PYTHIA code 

Usually, it is possible to fit hadron spectra up to high 
momentum, say pr ~ 2-3 GeV/c, within hydrodynamics 
by adjusting kinetic freeze-out temperature T th which is 
a free parameter in the model |3Jj. Thus it is unclear 
which value of T th should be used when one wants to 
add jet components into hydrodynamic components for 
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the description of high px part. However, we are free 
from this problem thanks to inclusion of the early chem- 
ical freeze-out picture into hydrodynamics. One of the 
authors studied the effects of chemical freeze-out temper- 
ature T ch which is separate from kinetic one T th in hy- 
drodynamic model in Ref. [2<| . It was found that the pr 
slope for pions remains invariant under the variation of 
T th and that the hydrodynamic model with early chem- 
ical freeze out is able to fit the transverse momentum 
distribution of pions up to 1-2 GeV/c. Therefore, it is 
certain to incorporate hard partons into the hydrody- 
namics with early freeze out in order to account for the 
high transverse momentum part of the hadronic spec- 
trum. We note that, since we do not assume thcrmaliza- 
tion for the high px jets, a hydrodynamical calculation 
with the initial conditions taken from pQCD+fmal state 
saturation model |37| is different from ours. 

In this paper, we shall study identified hadron spectra 
from low to high px within the hydro+jet model. In par- 
ticular, we focus on the influence of the hydrodynamic 
radial flow on the pQCD predictions for the transverse 
spectra. Parameters in the hydrodynamic part of the 
model have been already fixed by fitting the pscudorapid- 
ity distribution. Parameters related to the propagation 
of partons are also obtained by fitting the neutral pion 
suppression factor by PHENIX and are found to be con- 
sistent I38l with the back-to-back correlation data from 
STAR lU . 

The paper is organized as follows. In Sec. [n] we de- 
scribe the main features of our model. We will represent 
results of transverse momentum distributions for pions, 
kaons, and protons in Sec. IIII Al Nuclear modification 
factor (suppression factor) for identified hadrons and par- 
ticle ratios are discussed in Sec IIII Bl Elliptic flow for 
identified hadrons is discussed in Sec IIII Cl Section llVl 
summarizes this paper. 



II. MODEL DESCRIPTION 

In this section, we explain in some detail the hydro+jet 
model as a dynamical model to describe rclativistic heavy 
ion collisions. 



A. Hydrodynamics 

Let us start with the review of our hydrodynamics. 
Main features of the hydrodynamic part in the hydro+jet 
model are the following. 

Although initial conditions and pre-thermalization 
stages are very important sub ject s in the physics of heavy 
ion collisions (see, e.g., Ref. |40l l4l|). these are beyond 
the scope of this paper. Instead, assuming local ther- 
mal equilibrium of partonic/hadronic matter at an ini- 
tial time ro, we describe afterward the space-time evolu- 
tion of thermalized matter by solving the equations for 



energy-momentum conservation 

dyT^" = 0, = (e + P)u^u v - Pcf (1) 

in the full three-dimensional Bjorken coordinate 
(t, x,y,r) a ). Here e, P, and u M are, respectively, energy 
density, pressure, and local four velocity, r = \Jt 2 — z 2 
is the proper time and ?7 S = (1/2) ln[(i + z)/(t — z)] is the 
space-time rapidity. Throughout this paper, we consider 
baryon free matter ?ib = at RHIC energies. In order 
to obtain reliable solutions of Eq. especially in the 
longitudinal direction at collider energies, r and % are 
substantial choices for time and longitudinal directions 
rather than the Cartesian coordinate. 

Assuming N{ = 3 massless partonic gas for the QGP 
phase, an ideal gas EOS with a bag constant = 247 
MeV is used in the high temperature phase. We use 
a hadronic resonance gas model with all hadrons up 
to A(1232) for later stages of collisions. Possible fi- 
nite baryonic effects such as a repulsive mean field [42l | 
are not included because of the low baryon density at 
RHIC 43]. Phase transition temperature is set to be 
T c = 170 MeV. For the hadronic phase, a partial chem- 
ical equilibrium (PCE) model with chemical freeze-out 
temperature T ch = 170 MeV is employed to describe the 
early chemical freeze-out picture of hadronic matter. Al- 
though chemical freeze-out temperature T ch (~ 160-170 
MeV) is usually found to be larger than kinetic freeze-out 
temperature T th (^ 100-140 MeV) from statistical model 
analyses and thermal model fitting |44|, the sequential 
freeze out is not considered so far in the conventional hy- 
drodynamics except for a few work [2^, EE D3 ■ As a 
consequence of this improvement, the hadron phase cools 
down more rapidly than the one in usual hydrodynamic 
calculations in which T ch = T th is assumed [29ll45j . It 
should be emphasized that the slope of pions in the trans- 
verse momentum distribution becomes insensitive to the 
choice of the kinetic freeze-out temperature T th and that 
the hydrodynamics with early chemical freeze out repro- 
duces the RHIC data of the pion transverse momentum 
only up to 1.5 GeV/c 0- This is one of the strong mo- 
tivations which leads us to combine our hydrodynamics 
with non-thermalized hard components. 

From hydrodynamic simulations, we evaluate hadronic 
spectra which originate from thermalized hadronic mat- 
ter. For hadrons directly emitted from freeze-out hyper- 
surface S, we calculate spectra through the Cooper-Frye 
formula H3 

E dNi _ di f p^dcr^ 
~Wp~ ~ (2^r)3 h exp[(p^ - ^)/T th ] + 1 ' 

where di is a degeneracy factor, Hi is a chemical potential, 
p M is a four momentum in the center of mass frame of 
colliding two nuclei, and — (+) sign is taken for bosons 
(fermions). We should note the existence of chemical 
potentials [ii for all hadrons under consideration due to 
early chemical freeze out. Typical values at T th = 100 
MeV are as follows: fx n = 83 MeV, [i K = 181 MeV, 
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and Hp — /ip — 349 McV. For hadrons from resonance 
decays, we use Eq. (0) for resonance particles at freeze 
out and afterward take account of decay kinematics. Here 
these resonances also have their own chemical potentials 
at freeze out. We call the sum of the above spectra the 
soft component or the hydro component throughout this 
paper. 

Initial energy density at To = 0.6 fm/c is assumed to 
be factorized 



e{x,y,r) s ;b) 



tW(x,y;b)H( Vs ). 



(3) 



Here the transverse profile W(x, y; b) is proportional 
to the number of binary collisions and normalized as 
VF(0, 0; 0) = 1, whereas longitudinal profile H{j) s ) is flat 
and unity near midrapidity and falls off smoothly at large 
rapidity. In H{rj s ) 1 we have two adjustable parameters 
ryflat and 77Gauss which parametrize the length of flat re- 
gion near midrapidity and the width of Gaussian in the 
forward/backward rapidity region, respectively. These 
parameters are chosen so as to reproduce the shape of 
dN/dr) or dN/dY. 

We choose e max = 40 GeV/fm 3 , ?7fl at = 4.0, and 
?7Gauss = 0.8. As shown in Fig. ^ the pseudorapidity 
distribution of charged hadrons in 5% central collisions 
observed by the BRAHMS Collaboration j^j is satisfac- 
tory reproduced by using the above parameters. Here we 
choose an impact parameter as b = 2 fm for this central- 
ity These initial parameters give us an average initial 
energy density about 5 GeV/fm 3 in the transverse plane 
i] s = at r = 1 fm/c |5l| . A contribution from minijets 
is neglected in the hydrodynamic fitting, since it is less 
than 5% effect to the total hadron yield at RHIC when we 
define minijets as particles with transverse momentum 
larger than 2 GeV/c. Initial conditions for transverse 
profile are scaled by the number of binary collisions. It is 
found that the 20-30% semicentral collision data is also 
reproduced simply by choosing b as 7.2 fm in the trans- 
verse profile W |52j . 

In Fig. [3 wc show the transverse spectra for negative 
pions, negative kaons, and protons in Au+Au collisions 
at ^sjvtv = 200 GeV from the hydrodynamic model for 
impact parameters b = 2.0 fm and 7.2 fm. Thermal 
freeze-out temperature T th = 100 MeV is used in the 
calculation. This choice is consistent with the data at 
^snn = 130 GeV [2^. The flatter behavior at lowpr for 
kaons and protons is indeed a consequence of the radial 
flow effect. A remarkable feature on the hydrodynamical 
result is that p/n~ > 1 and K~ /tt~ ~ 1 above pr ~ 2 
GeV/c. It is, however, questionable to assume thermal- 
ization at high p^ region. In fact, hydrodynamical pre- 
dictions overestimate elliptic flow data at the large trans- 
verse momentum region. It is interesting to ask at which 
Pt hydrodynamic behavior ceases and switches to pQCD 
results. We will see in the next section how these hydro- 
dynamical results are modified by including the pQCD 
hard component. 
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FIG. 1: (Color online) Pseudorapidity distribution of charged 
particles in Au + Au collisions at y / sjvjv = 200 GeV is com- 
pared to data from BRAHMS l50l . Solid (dashed) line repre- 
sents the hydrodynamic result at b=2.0 (7.2) fm. 
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FIG. 2: (Color online) Transverse momentum spectra for 
negative pions, negative kaons, and protons from the hydro 
model with early chemical freeze out in Au+Au collisions at 
y^ijvjv = 200 GeV. We choose an impact parameter as b — 2.0 
(7.2) fm corresponding to 0-5% (20-30%) centrality. Yields 
are divided by 10 3 for b — 7.2 fm results. 



B. Jet propagations 

For the hard part of the model, we generate hard par- 
tons according to a pQCD parton model. The number of 
jets at an impact parameter b are calculated from 



NhardW 



d 2 r ± a jct T A (r ± - b/2)T B (r ± + 6/2), (4) 



where a )C t is a hard cross section from leading order 
pQCD convoluted by the parton distribution functions 
and multiplied by a iGfactor which takes into account 
higher order contributions. Ta and r± are, respectively, a 
nuclear thickness function normalized to be J d 2 r±T A = 
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A and a transverse coordinate vector. Here we use the 
Woods-Saxon distribution for the nuclear density profile. 
We use PYTHIA 6.2 for the generation of momen- 
tum spectrum of jets through 2^2 QCD hard processes. 
Initial and final state radiations are used to take into 
account the enhancement of higher-order contributions 
associated with multiple small-angle parton emission. 

Scale Q 2 dependent nuclear shadowing effect is in- 
cluded for the mass number A nucleus assuming the im- 
pact parameter dependence |5^ |: 



S{A, x, Q 2 ,r ± ) = l + [S(A, x, Q 2 ) - 1]- 



AT A ( ri _) 



l fd 2 r ± T A (r ± ) 2 ' 
l-l (5) 

where the EKS98 parametrization [54j is used for 
S(A, x, Q 2 )- Then the nuclear parton distribution func- 
tion in this model has the form 

/ A (A,x,Q 2 ,r ± ) = S(A, Xl Q 2 ,r ± ) 

jf P (x,Q 2 )+ ( ^-J^-fn(x,Q 2 ) 

(6) 

where f p (x, Q 2 ) and f n (%, Q 2 ) are the parton distribution 
functions for protons and neutrons. We simply assume 
the charge of a nucleus to be Z = A/2 in consistency with 
the soft part, since our fluids are assumed to be isospin 
symmetric as well as baryon free matter. 

Cronin effect [5j|, which has also been discovered in 
recent RHIC experiments [5(|, is usually considered as 
the multiple initial state scattering effect. Understand- 
ing this effect becomes an important subject in RHIC 
physics |52i E3, ■ We employ the model in Ref. J5?J to 
take into account the multiple initial state scatterings, in 
which initial fc<r is broadened proportional to the number 
of scatterings: 

(k 2 T ) NA = (k 2 T ) NN + 5 2 {Q 2 ) (a NN T A (r x ) - 1) , (7) 

where <7nn is the inelastic nucleon- nucleoli cross sec- 
tion and 5 2 (Q 2 ) is the scale dependent kr broadening 
per nucleon-nucleon collision whose explicit form can be 
found in Ref. [57j . 

We need to specify a scale which separates a soft sec- 
tor from a hard sector, in other words, a thermalized 
part from a non-thermalized part in our model. We in- 
clude minijets with transverse momentum prjet larger 
than 2 GeV/c just after hard scatterings in the simula- 
tion. These minijets explicitly propagate through fluid 
elements. 

Since we only pick up high pr partons from PYTHIA 
and throw them into fluids, there is ambiguity to con- 
nect color flow among partons. Thus we use an inde- 
pendent fragmentation model option in PYTHIA to con- 
vert hard parton to hadrons instead of using the default 
Lund string fragmentation model. We note that the in- 
dependent fragmentation model should not be applied at 
low transverse momentum region. We have checked that 
the neutral pion transverse spectrum in pp collisions at 



RHIC : 6Q| is well reproduced by selecting the if-factor 
K = 2.5, the scale Q = prjct/2 in the CTEQ5 lead- 
ing order parton distribution function [6]]], and the pri- 
mordial transverse momentum (fc T )jv\2V = 1 GeV 2 /c 2 as 
shown in Fig [JJ As shown in the bottom panel of the 
Fig. [31 independent fragmentation model predictions for 
pions and kaons are very close to those from the Lund 
string fragmentation model in pr > 2 GcV/c, where 
K = 2 and Q = pxjct/2 is used in the Lund string model 
case and non-perturbative inelastic soft processes are in- 
cluded. However, the yield of protons from the indepen- 
dent fragmentation scheme becomes much less than that 
from Lund string model predictions as seen in Fig. |3J We 
found that the Lund fragmentation scheme is favored in 
terms of the recent STAR data of protons in pp colli- 
sions |62l |. In what follows, we make corrections for our 
pr spectra of kaons and protons in Au+Au collisions ac- 
cording to the result in the bottom panel of Fig. [3J 

In order to see the theoretical uncertainties on the frag- 
mentation scheme deeply, we also plot the results from 
NLOpQCD calculations with the MRST99 H set 
of parton distribution functions. In Fig. UJ we show re- 
sults from two different fragmentation functions. The 
solid lines are obtained from KKP fragmentation func- 
tions |65| with renormalization scale fj>, factorization scale 
M, and fragmentation scale Mf equal to pr- NLOpQCD 
prediction with KKP fragmentation functions is consis- 
tent with the pion data. NLOpQCD predictions with the 
Kretzer fragmentation functions j6|| assuming /i = M = 
Mf = pr/2 underestimate pion yields, while yields for 
kaons and protons are the same as the predictions from 
the PYTHIA default Lund string fragmentation model. 

Initial transverse positions of jets at an impact param- 
eter b arc determined randomly according to the proba- 
bility P(r±, b) specified by the number of binary collision 
distribution, 



P(rx, 6) <x T A (rj_ + b/2)T A ( r± _ - 6/2). 



(8) 



Initial longitudinal position of a parton is approximated 
by the boost invariant distribution |67j : % = Y , where 
Y = (l/2)ln[(E + p z )/(E ~ p z )] is the rapidity of a par- 
ton. Jets are freely propagated up to the initial time To 
of hydrodynamic simulations by neglecting the possible 
interactions in the pre-thermalization stages. Jets are 
assumed to travel with straight line trajectory in a time 
step: 



An = 



Pi 



niT cosh(Y — i] s ) 



1 



A?y s = — tanh(Y — 7y s )Ar, 
r 



(i = x,y), (9) 
(10) 



where wit = y m 2 + p 2 -, is a transverse mass. 

Jets can suffer interaction with fluids and lose their 
energies. We employ the approximate first order for- 
mula (GLV formula) in opacity expansion from the reac- 
tion operator approach Q for the energy loss of partons 
throughout this work. The opacity expansion is relevant 
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FIG. 3: (Color online) Comparison with various models for 
inclusive pion, kaon, and proton transverse momentum dis- 
tributions in pp collisions at a/s = 200 GeV. Solid and dot- 
ted histograms correspond to the results from PYTHIA with 
independent fragmentation and default Lund fragmentation 
respectively. Solid and dotted lines are, respectively, from 
NLOpQCD calculations with KKP and Kretzer fragmenta- 
tion functions. 



for the realistic heavy ion reactions where the number of 
jet scatterings is small. The energy loss formula for co- 
herent scatterings in matter has been applied to analyses 
of heavy ion reactions taking into account the expansion 
of the system |32, E^. The approximate first 

order formula in this approach can be written as 



AE = C I drp (r, x (r)) (t - r ) In 
'to 



2E 
p?L 



(11) 



Here C is an adjustable parameter and p(r, x) is a ther- 
malized parton density in the local rest frame of fluid el- 
ements in the hydro+jet approach (6^|. x (r) and Eq are 
the position and the initial energy of a jet, respectively. 
The initial energy Eq in Eq. (jllfl is Lorentz-boosted by 
the flow velocity and replaced by PqU^ where pg is the 
initial four momentum of a jet and is a local fluid ve- 
locity. We take a typical screening scale p = 0.5 GeV and 
effective path length L = 3 fm which is chosen from the 
lifetime of the QGP phase. Here we choose C = 0.45 [g^ 
which is found to reproduce the neutral pion Raa defined 
by Eq. l|12f) ■ Our purpose here is not a detailed study 
of jet quenching mechanisms. Instead, we first fit the 
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FIG. 4: (Color online) Jet quenching rate Nj ct (T) / Nj et (To) for 
Pt — 5 GeV/c jets in Au+Au collisions at ^/sNN = 200 GeV. 
Jet quenching rate for 10 GeV/c jets is very similar to that 
of the 5 GeV/c jets. 



suppression factor for neutral pions and next see other 
hadronic spectra. 

Feedback of the energy to fluid elements in central col- 
lisions was found to be about 2% of the total fluid energy. 
Hence we can safely neglect its effect on hydrodynamic 
evolution in the case of the appropriate amount of energy 
loss. 

In Fig. 01 we show the jet quenching rate as a function 
of proper time for 5 GeV/c jets. We count the num- 
ber of partons with 4.5 < PTjct < 5.5 GeV/c at each 
time step, and then define the ratio of the current num- 
ber of jets to the initial number of jets Nj et (T) / Nj et (To) ■ 
Most jet quenching is completed at early times less than 
4 fm/c. For comparison, we also plot the jet quenching 
rate for a constant energy loss case dE/dx oc p{r). Jet 
quenching is almost finished at r ~ 2 fm/c in the case of 
constant energy loss. From Fig.|3| the degree of decrease 
for the jet quenching rate in the GLV formula becomes 
milder and continues longer than that in the incoherent 
model. This is due to the existence of t in the integrand 
in Eq. (|11() which comes from the property of coherent 
(Landau-Pomeranchuk-Migdal_j7ij|) effect. Contrary to 
the simple Bjorkcn's ansatz |67j. p(r) = PqTo/t, there 
exists transverse flow and the parton density profile in 
the transverse plane is not flat in our simulations. This 
is the reason why jets are quenched only in the QGP 
phase and why jet quenching in the mixed phase is to- 
tally negligible. 

We include p± broadening accompanied by the en- 
ergy loss of jets with the formula ~ / drp{r) as in 
Ref. jj^l ■ We found that this effect is small in all results 
in this paper. 

Within our model, we neglect energy loss before ther- 
malization, in our case, r < 0.6 fm/c. One would ask 
if it is important to take into account the energy loss 
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FIG. 5: (Color online) Transverse momentum spectra for neg- 
ative pions, negative kaons, and protons from the hydro+jet 
model in Au+Au collisions at ^/snn = 200 GeV at the im- 
pact parameter of b = 2.0 and 7.2 fm. Yields are divided by 
10 3 for 6 = 7.2 fm results. 



effects before thermalization because parton density has 
the maximum value. We can, however, fit the suppres- 
sion factor Raa by rescaling the energy loss parameter 
C when the initial time To is changed. The question 
about the jet quenching before thermalization is beyond 
our model description. As a possible model for a study 
of jet interactions at early times, propagation of jets in 
the classical Yang-Mills fields based on the idea of the 
Color Glass Condensate 0, IxU is proposed in Ref . [f2 • 
It would be interesting to take a numerical results from 
the full lattice calculations for the calculations of jet 
energy loss at the very early stages of the collisions. 



III. RESULTS 

We discuss in this section transverse dynamics for pi- 
ons, kaons, and protons from the hydro+jet model fo- 
cusing on the intermediate px where interplay between 
soft and hard components is expected to be crucial. As 
mentioned in the previous section, a parameter for jet 
quenching C was already fixed by fitting the observed 
data for neutral pions in central Au+Au collisions from 
PHENIX. Freeze-out temperature T th = 100 MeV is used 
for hydrodynamics. All results in this section are for 
midrapidity | r\ |< 0.35. 



A. Transverse momentum distributions for 
identified particles 

First, we show the transverse momentum distributions 
for pions, kaons, and protons from the hydro+jet model 
in Fig. [3] in central as well as semicentral Au+Au col- 
lisions at RHIC. Each spectrum is the sum of the soft 



component and the hard component. Before summation, 
the hard component is multiplied by a "switch" func- 
tion [3l| {1 + tanh[2(pT — Pr,cut)]}/2 (where px is in the 
unit of GeV/c and px,cut = 2 GeV/c) in order to cut the 
unreliable components from the independent fragmenta- 
tion scheme and also to fit Raa for neutral pions [jj. 
We have checked the cutoff parameter dependence in the 
switch function on the pion spectrum and found that 
we are not able to fit the pion data anymore even with 
PT.cut = 1-8 or 2.2 GeV/c. So the ambiguity of the cutoff 
can be removed to fit the pion data within our approach. 

At low transverse momentum region px < 1 GeV/c, 
the shapes remain the same as hydro predictions as one 
can check from Fig. [3 Also at high transverse momen- 
tum, spectra are identical to those of pQCD predictions 
with an appropriate amount of jet quenching. 

Our calculation includes interactions of minijets with 
QGP fluids. We also note that there remains a pQCD 
like power law behavior in all hadrons at high transverse 
momentum. This may indicate no hint for the thermal- 
ization at high transverse momentum. However, energy 
loss results in a parallel shift of hadronic spectra, since 
the energy loss model used in this paper shows almost flat 
quenching pattern as shown in our previous analysis |38l | . 

In Fig. we decompose the spectra into hydro parts 
and minijet parts. Here the yields from hard components 
are multiplied by the switch function again. It is seen 
that both soft and hard components are important for 
the hadron spectra in the transverse momentum of the 
range around 2<,px^S 5 GeV/c depending on the hadron 
mass. We can define the crossing point of transverse 
momentum px, cross at which the yield from the soft part 
is identical to that from the hard part, px, cross moves 
toward high momentum with mass of particles because of 
the effects of radial flow. In central collisions, px. cross ~ 
1.8, 2.5, and 3.5 GeV/c for pions, kaons, and protons, 
respectively. Minijet spectra are recovered at px ~ 3.4 
GeV/c for pions, px ~ 4.0 GeV/c for kaons, and px ~ 5.0 
GeV/c for protons. 

We give some remarks here: 

(a) The point at which hydrodynamic and pQCD spec- 
tra cross is determined by the dynamics of the system: 
The radial flow pushes the soft components toward high 
px region, while the dense matter reduces the pQCD 
components through parton energy loss. The crossing 
of two spectra causes by interplay of these two effects. 

(b) At px = 2-3 GeV/c, the yields of pions and kaons 
are no longer occupied by soft hydrodynamic component. 
On the other hand, the proton yield from pQCD predic- 
tion is about ten times smaller than that of hydro in the 
transverse momentum region. 

(c) One may try to extract the strength of radial flow 
and the kinetic freeze-out temperature from experimen- 
tal data through the hydrodynamics-motivated fitting 
model. Then one should pay attention to the fitting range 
of the transverse momentum. In particular, px spectrum 
for pions may have no room to fit by a simple thermal 
spectrum: Contribution from resonance decays becomes 
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FIG. 6: (Color online) Each contribution from hydrodynam- 
ics and minijets for n~ , K~ , and p in Au+Au collisions at 
^/sNN = 200 GeV at the impact parameter of b = 2.0 fm. 
Yield of negative kaons (protons) is divided by 10 3 (10 6 ). 
PHENIX data are from Ref. Fl. 
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FIG. 7: 



(Color online) Impact parameter dependence of the 
suppression factors Raa in Au+Au collisions at ^Jsnn = 200 
GeV as a function of pt for tt~ , K~ , and p. Raa for charged 
hadrons is also shown in dotted lines. Experimental data of 
Raa for neutral pions is obtained by PHENIX. For details, 
see text. 



models. 



important below pt ~ 0.5 GeV/c, while the hard compo- 
nent slides in the soft component near px ~ 1.0 GeV/c. 

(d) We predict positions of the inflection point where 
Pt spectrum becomes convex to concave: pt ~ 3 GeV/c 
for kaons and ~ 4 GeV/c for protons. These are the in- 
dicators of a transition from soft physics to hard physics. 

The amount of the hydrodynamic contributions to the 
hadron yields for each particle found in the hydro+jet 
model is very similar to that found in Ref. p5j in which 
hybrid parametrization of hydrodynamics with the spec- 
tral shape in pp collisions. It is also remarkable that 
bar yon junction [Til ITsj and quark coalescence mod- 
els |2Ct I2TL |22| predicts the same behavior. Quark coa- 
lescence models are successful in explaining the mass de- 
pendence of pt slopes |75ll76j . For example, one can eas- 
ily understand the difference of the transverse slopes of 
baryons and mesons from a quark coalescence hadroniza- 
tion mechanism. A baryon momentum is a sum of three 
quarks (quark momenta must be almost parallel in or- 
der to cluster), but a momentum of mesons is a sum of 
two quarks. It is interesting to see, for example, 4> me ~ 
son spectrum in order to distinguish the mass effects in 
hydrodynamics from meson-baryon effects in coalescence 



B. Suppression factors and particle ratios 

We now turn to the study of the suppression factors 
Raa for each hadron defined by 

d.N A+A 



It is very instructive to study Raa behaviors for iden- 
tified hadrons toward a comprehensive understanding of 
intermediate transverse momentum region. 

Figure [7| shows the suppression factors Raa for pi- 
ons, kaons, protons, and charged hadrons respectively 
in Au+Au collisions at RHIC for impact parameters 
b = 2.0, 5.5, and 7.2 fm. Our results for pions arc 
compared with PHENIX data @. We use pp spectra 
from Lund string model for protons and kaons in the 
plots. Raa for protons using the independent fragmen- 
tation model becomes too large ~ 2.5 at px ~ 2 GeV/c. 
This result simply comes from the fact that the indepen- 
dent fragmentation model is inconsistent with pp data 
for protons as discussed in Sec. [H] Note that the nu- 
merator in Eq. I|12|) is almost free from the hard compo- 
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FIG. 8: (Color online) Ratios of N p to N n - and N K - to N n - 
as a function of pt in Au+Au collisions at impact parameter 
6 = 2 fm. PHENIX data [74j are also ploted for comparsion. 



nents at pt < 3 GeV/c in proton case. We find protons 
are not suppressed Raa > 1 at a momentum range of 
1.5 < pt < 2.5 GeV/c. Pions, on the contrary, are largely 
suppressed for all momentum range. Our calculations for 
protons become identical to those of pQCD predictions 
at a momentum above 5 GeV/c. This is the same re- 
sult as other model predictions 0,H2,|2f|. In any case, 
these results are easily understood from Fig.[S] The cross- 
ing point pt, cross depends on the hadronic species, thus 
Raa only for pions reflects jet quenching effect, while 
the larger value of Raa for protons simply comes from 
radial flow, not the absence of jet quenching. We should 
mention that above pt ~ 5 GeV/c, suppression factors 
for identified hadrons converge to almost the same value. 
It is also seen that the suppression factors for kaons and 
protons have almost no ccntrality dependence within this 
impact parameter range. 

Recent data from PHENIX P3 and STAR Q3 for pro- 
tons and A's show that the nuclear modification factors 
for p, p, and A in the pr range of 1.5 < pr < 4.5 GeV/c 
are almost constant. However, our results of Raa for 
protons seem to decrease to smaller value with transverse 
momentum faster than data. 

From RHIC data USQUE!, Raa for char S ed P arti " 
cles is larger than the one for pions in moderate high pt 
region. In our model, this also simply results from the 
average of the above three suppression factors weighted 
by each yield (see also Fig. a s shown in Fig. [7] by the 
dotted lines. 

We show in Fig. [8] proton to negative pion ratio and 
negative kaon to negative pion ratios as a function of the 
transverse momentum in Au+Au collisions at ^/sJjn = 
200 GeV for the impact parameter of b = 2 fm together 
with the PHENIX data |z3 - Without depending on 
baryon junction mechanism or quark coalescence models, 
we also obtain that p/^~ ratio becomes close to unity due 



to the consequences of hadron species dependent pt, cross- 
Ratios become p/tt~ ~ 0.2 and K~ /ir~ ~ 0.3 above 
Pt ~ 5 GeV/c which are the consequences of pQCD pre- 
dictions. It should be noted that, if the baryonic and 
isospin chemical potentials are included in the hydro- 
dynamic simulation, p/ir~ ratio can slightly be changed 
in low pt region: Baryon (isospin) chemical potential 
pushes up (down) proton yield from hydrodynamic com- 
ponents. 



C. Elliptic flow for identified particles 

Azimuthal asymmetry for non-central heavy ion colli- 
sions is generally considered to be generated only by the 
final state interactions of matter created in the collisions. 
In hydrodynamic models, elliptic flow is created by the 
anisotropic initial configuration of high pressure matter 
which might be the QGP phase. 

Hydrodynamic predictions on the transverse momen- 
tum dependence of elliptic flow V2 show almost linear in- 
crease for all particles. However, the experimental data 
saturate at highpT [3,|23,lz3- More interestingly, pion 
i>2 is larger than that of protons at pt < 1 GeV/c, while 
proton V2 becomes larger than pion V2 at some pt |79j| . 
Hydrodynamic calculations are successful in reproduc- 
ing the mass dependence of the V2 in the low transverse 
momentum region j2|| IH, |2jj . w 2 for pions are al- 
ways greater than that of protons in hydrodynamics and, 
eventually, i>2 becomes almost mass independent at high 
transverse momenta as shown in Fig. On the other 
hand, to understand the observed azimuthal asymme- 
try at large pt, it was showed that the jet interaction 
with matter generates the azimuthal asymmetry for non- 
central collisions [M HI 113 

We demonstrate in Fig. [5] that, by combining mini- 
jet components with hydrodynamics, pion V2 can be re- 
duced faster than proton i>2 at moderate high transverse 
momentum. The hydro+jet predictions on V2 for iden- 
tified particles in Au+Au collisions at RHIC for impact 
parameter b = 5.5 fm are compared to hydro results in 
Fig. [5J The magnitude of «2 for kaons and protons be- 
comes larger than u 2 for pions at about pt > 1.3 GeV/c. 
The shape of i>2 for pions saturates faster than those of 
kaons and protons, because the fraction of hydro compo- 
nents for pions are much smaller than that for kaons or 
protons in this pt region. This is again the consequence 
of radial flow effect. We demonstrate that the saturation 
point in transverse momentum depends on the hadron 
mass. As a whole effect of the sum of pions, kaons, and 
protons, the saturation point of V2 for charged particles 
in transverse momentum is turned out to be pt = 1.5 
GeV/c in our model at b = 5.5 fm as one can read from 
Fig. ED 

Our semi-macroscopic model produces consistent 
behavior in V2 with the experimental data from 
PHENIX [79|. Recently, a microscopic description of 
quark coalescence model shows the crossing of mc- 
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va ~ v<j>- Therefore, it is interesting to see, for example, 
</> meson elliptic flow to clarify the origin of the elliptic 
flow. 

We have studied vi in the momentum range where both 
soft and hard contributions are important. It is interest- 
ing to see V2 up to 10 GeV/c. Experimental data show 
that vi at high momentum saturates |77j . Systematic 
study on the elliptic flow within our model is under way 
including centrality as well as rapidity dependence. 
IV. SUMMARY 

We have studied the interplay of soft and hard com- 
ponents by looking at pt spectra, suppression factors, 
hadron ratios, and elliptic flow for identified particles 
within the hydro+jet model. By taking into account both 
hydrodynamic radial flow and quenched pQCD spectra, 
it was found that pr, cross, at which the yield from the soft 
component is identical to the one from the hard compo- 
nent, depends on the hadron species: pr. cross ~ 1.8,2.5, 
and 3.5 GeV/c for pions, kaons, and protons in Au+Au 
central collisions at RHIC. This difference comes from 
the interplay between the radial flow for the soft part 
and the jet quenching for the hard part. From the conse- 
quences of the interplay between soft and hard hadronic 
components, we showed p/-n~ ~ 1 and Raa(pt) > 1 at 
intermediate pr for protons. We also showed that the 
mass dependence of the strength of i>2(pr) in the inter- 
mediate pt region is also explained by the radial flow 
+ pQCD components. Hydrodynamic radial flow plays 
an important role to understand the transverse dynamics 
when hadron mass is large. 



FIG. 9: (Color online) Fransverse momentum dependence of 
the elliptic flow V2 (pt) for 7r _ , K~ , and p in Au+Au collisions 
at y/SNN = 200 GeV at impact parameter b — 5.5 fm from the 
hydrodynamic model (upper panel) and the hydro+jet model 
(lower panel). V2 for charged hadrons are also represented in 
dotted lines. 

son and baryon t^'s at pr ~ 1 GeV/c. In the simple 
coalescence model where all partons have similar elliptic 
flow, elliptic flow for baryons roughly 1.5 times stronger 
than for mesons. On the other hand, our approach will 
have only mass dependence on the elliptic flow indicating 
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